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Abstract 

We show that a TeV scale inverse seesaw model for neutrino masses can be 
0^ ■ realized within the framework of a supersymmetric 50(10) model consistent with 

gauge coupling unification and observed neutrino masses and mixing. We present 
our expectations for non-unitarity effects in the leptonic mixing matrix some of 
which are observable at future neutrino factories as well as the next generation 
searches for lepton flavor violating processes such as \x — > e + 7. The model has TeV 
scale Wr and Z' bosons which are accessible at the Large Hadron Collider. 

1 Introduction 

A precise understanding of the origin of observed neutrino masses and mixing is one of 
the major goals of particle physics right now. A simple paradigm for understanding the 
smallness of the masses is the seesaw mechanism [1] where one introduces three Standard 
Model (SM) singlet right-handed (RH) neutrinos with Majorana masses M^, which mix 
with left-handed (LH) ones via the Yukawa coupling LHN . The resulting formula for 
light neutrino masses is given by Ai u = —M^M^M^, where is the Dirac neutrino 
mass. Since the SM does not restrict the Majorana mass M N , we could choose this to 
be much larger than the weak scale thereby providing a natural way to understand the 
tiny neutrino masses. This is called the type I seesaw. There are several variations of this 
mechanism where one replaces the RH neutrino by either a SM triplet Higgs field (type II 
seesaw) [2] or SM triplet of fermions (called type III seesaw) [3]. A great deal of attention 
has been devoted to testing these ideas. As far as the type I seesaw is concerned, the 
prospects of testing this depends on the scale M N as well as any associated physics that 
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comes with it at that scale. It can be accessible to current and future collider experiments 
if the scale is not far above a TeV. A different way to test the type I seesaw mechanism 
follows from the observation that this mechanism involves the mixing of the LH neutrinos 
with SM singlet heavy neutrinos as a result of which there would in general be violation of 
unitarity of the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing matrix that describes 
only the mixing of the three light neutrinos. One could contemplate searching for such 
effects in oscillation experiments pEj. However, in the type I seesaw case, the resulting 
mixing effects are of order j^- and since neutrino masses are in the sub-eV range, such 
non-unitarity effects are too small to be observable for generic high scale seesaw models. 
This would also be true with TeV mass RH neutrinos unless there are cancellations to get 
small neutrino masses from large Dirac masses using symmetries (see for example cases 
with [5J). We note that the non-unitarity effects are also there in the type III seesaw case 
but not in the type II case even though they have other interesting effects such as LFV 
processes [6]. 

Since the testability of seesaw is intimately related to the magnitude of the 
seesaw scale, a key question of interest is whether there could be any theoretical guidelines 
for the seesaw scale. In such a case, the searches for seesaw effects in experiments could 
then be used to test the nature of physics beyond the standard model. It is well-known 
that [7] the simplest grand unified theory (GUT) realizations of the seesaw mechanism are 
based on the SO (10) group which automatically predicts the existence of the RH neutrinos 
(along with the SM fermions) required by the seesaw mechanism. An advantage of GUT 
embedding of the seesaw mechanism is that the constraints of GUT symmetry tends 
to relate the Dirac neutrino mass to the charged fermion masses thereby making a 
prediction for the seesaw scale Mjv from observations. For type I seesaw GUT embedding, 
typical values for the are very high (in the range of 10 10 - 10 14 GeV). This makes 
both the collider as well as non-unitarity probes of seesaw impossible. The key feature 
that leads to such restrictions in type I seesaw case is the close link between the B — L 
breaking RH neutrino mass and the smallness of the LH neutrino masses. 

A completely different realization of the seesaw mechanism is the so-called 
inverse seesaw mechanism j8], where instead of one set of three SM singlet fermions, 
one introduces two sets of them N iy Si (i = 1,2,3). In the context of 5*0(10) models, 
since one of the two sets can be identified with the SM singlet neutrino in the 5*0(10) 
16-representation containing matter, the others would have to be a separate set of three 
50(10) singlet fermions. Due to the existence of the second set of singlet fermions (and 
perhaps additional gauge symmetries e.g. 50(10)), the neutrino mass formula in these 
models has the form 

m u ~ M D M^(i (Ml) ~ l M T D = FfiF T (1) 

where \x breaks the lepton number. Because of the presence of this new mass scale in this 
theory, the seesaw scale can be very close to a TeV even for "large" Dirac masses. This 
makes the tests of this possibility in colliders much more feasible. In fact it has recently 
been argued that [9] the inverse seesaw scenario can also lead to non-negligible non- 
unitarity effects which can be accessible at the future long-baseline neutrino oscillation 
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experiments. There are also significant lepton flavor violation (LFV) effects in these 
models as noted many years ago in Ref. [TO]. These possibilities have generated a great 
deal of interest in the inverse seesaw models in recent days [11]. Our effort in this paper 
focuses on possible grand unification of inverse seesaw models. Similar unification studies 
have been performed in Ref. [12], but they have not addressed the non-unitarity issues. 

An interesting question is whether such models are necessarily compatible with 
grand unification when the seesaw scale is in the TeV range and if so what kind of non- 
unitarity effects they predict. We find that it is indeed possible to embed the TeV scale 
inverse seesaw models within a simple SO(10) framework consistent with gauge coupling 
unification and realistic fermion masses. The 50(10) symmetry helps to reduce the num- 
ber of parameters in the inverse seesaw matrix, once we require degeneracy of the TeV 
scale RH neutrinos to have successful resonant leptogenesis. Within this set of assump- 
tions, we present our expectations for the non-unitarity effects as well as consequences for 
lepton flavor violation which are in the testable range in future experiments. 

This paper is organized as follows: In Section 2, we describe the general frame- 
work of the inverse seesaw model and its embedding into a generic supersymmetric .SO (10) 
theory. In Section 3, we analyze the non-unitarity predictions of the inverse seesaw model. 
In Section 4, we investigate a specific S'O(IO) breaking chain and obtain the gauge cou- 
pling unification with TeV scale left-right symmetry with unification scale consistent with 
proton decay bounds. In Section 5, we analyze the renormalization group (RG) evolution 
of the Yukawa couplings and obtain the running masses for quarks and leptons at the 
unification scale to check that our model leads to realistic fermion masses. In Section 6, 
we determine the Dirac neutrino mass matrix using the results of Section 5. In Section 7, 
we study the implication of our model on non-unitarity effects and its phenomenological 
consequences. A brief summary of the results is presented in Section 8. In Appendix A, 
we have given the expressions for the masses of the 5*0(10) Higgs multiplets in our model, 
and in Appendix B, we have derived the RG equations for the quark and lepton masses 
and the CKM mixing elements in the context of a supersymmetric left-right model. 

2 The Inverse Seesaw Model 

The inverse seesaw scheme was originally suggested [8] for theories which lack the repre- 
sentation required to implement the canonical seesaw, such as the superstring models. As 
noted in the introduction, the implementation of the inverse seesaw requires the addition 
of three extra SM gauge singlets Si coupled to the RH neutrinos iVj through the lepton 
number conserving couplings of the type NS while the traditional RH neutrino Majorana 
mass term is forbidden by extra symmetries. The lepton number is broken only by the 
self coupling term /iS" 2 . The mass part of the neutrino sector Lagrangian in the flavor 
basis is given by 

£ mass = (VM D N + NM N S + h.c.) + SfiS (2) 

where fi is a complex symmetric 3x3 matrix (with dimension of mass), and Mp and 
are generic 3x3 complex matrices representing the Dirac mass terms in the v - N and N 
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- 5 sectors respectively. In the basis {v, N, 5}, the 9x9 neutrino mass matrix becomes 




M v = Ml M N (3) 



The LH neutrinos can be made very light (sub-eV scale), as required by the oscillation 
data, even for a low M N , much smaller than the unification scale (M N <C M G ), provided 
/i is sufficiently small, \x -C M N , as the lepton number breaking scale \i is decoupled from 
the RH neutrino mass scale. Assuming fi <C <C Mjv (with Mn ~ TeV), the structure 
of the light neutrino Majorana mass term at the leading order in M^M^ 1 is given by 
Eq. (JTJ, where F = M^M^ 1 is a complex 3x3 matrix. We note that in the limit /i — > 
which corresponds to the unbroken lepton number, we have massless LH neutrinos as in 
the SM. In reality, a small non- vanishing \x can be viewed as a slight breaking of a global 
£7(1) symmetry; hence, the smallness of \x is natural, in the 't Hooft sense [13], even 
though there is no dynamical understanding of this smallness. 

The generic form of the inverse seesaw matrix in Eq. ([3]) has more parameters 
than the usual type I seesaw. However, if we embed this theory into a grand unified theory 
such as SO (10), that will help in reducing the parameters as we show below. In order to 
embed the inverse see-saw mechanism into a generic 5*0(10) theory, we have to break the 
B — L symmetry by using a 16 © 16 pair rather than the 126 © 126 pair of Higgs 
representation. All the SM fermions are accommodated in a single 16^ representation of 
50(10) and we use three copies of 16' F for three generations. For each of them, we add 
a gauge singlet fermion V F to play the role of 5j. We assume more than one copy of 10# 
Higgs multiplets in order to have a realistic fermionic spectrum. 

The 50(10) invariant renormalizable Yukawa superpotential is given by 

w Y = /^.i6* F i6 J F io^ + f^w^Wv + / /,/r / /i;, (4) 

After the B — L symmetry breaking, we get the neutrino mass matrix in Eq. ([3]) with 
Mo = hv u and M^ = fvR, where v u is the vacuum expectation value (VEV) of one (or, 
a linear combination) of the 10#'s and vr the VEV of the 16#. In a typical TeV-scale 
scenario with v u ~ 100 GeV (electroweak scale) and Vr ~ TeV, assuming \i <C v u < Vr, 
we find the lightest neutrino mass from Eq. ([1]) in a one generation theory to be 

m„ * „ g£) (5) 

and the two other heavy eigenstates with mass of order fvR. Thus, we can get sub-eV 
light neutrino mass for \i ~ keV. Since this is a super symmetric theory, such small values 
do not receive radiative corrections and keep the model natural. In the following section, 
we consider three generations which then results in the non-unitarity effect. 

It is important to note that in our model, we do not need to impose a dis- 
crete i?-parity to our matter fermions, unlike the usual 16# 5O(10) models discussed in 
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literature, in order to prevent fast proton decay via dimension-4 operators of the type 
-^16^16^16^16^ because these operators are already suppressed by a factor ~ 
10~ 15 for a low-scale B — L breaking with (16 h) ~ TeV. 



3 Non-unitarity Effects 

The 3x3 light neutrino mass matrix in Eq. (DD) can be diagonalized by a unitary trans- 
formation: 

U^raJJ* = m u = diag(mi, m 2 , m 3 ) (6) 

where U is the standard PMNS matrix. Since the above diagonalization of m v does 
not diagonalize the matrices and fi, there will be off-diagonal mixing between the 
different light neutrinos even after diagonalization of m u due to their mixing with the 
heavy neutrinos. In other words, in the basis where the charged-lepton mass matrix is 
diagonal, U is only a part of the full mixing matrix responsible for neutrino oscillations. 
We have to examine the full 9x9 unitary matrix V which diagonalizes the mass matrix 
M v given by Eq. (j3J): 

V^M U V* = M„ = diag(m h m Nj ,m^ k ) (i, j, k = 1, 2, 3) (7) 

We can decompose V into the blocks 

V = 

Then the upper-left sub-block V3X3 will represent the full (non-unitary) PMNS mixing 
matrix. For a TeV-scale Mjv and a reasonably small /1, it is sufficient to consider only up 
to the leading order in F. Then the new PMNS matrix becomes 




M= V 3 x3- (l- l -FF^U (9) 

In the commonly used parametrization [15J, M = (l—i])U, and hence, all the non-unitarity 
effects are determined by the Hermitian matrix r\ ~ \FF^ which depends only on the 
mass ratio F = M^M^ 1 and not on the parametrization of the PMNS matrix. 

The LH neutrinos entering the charged- current interactions of the SM now 
become superpositions of the nine mass eigenstates (i>i,Ni,Ni) and at the leading order 
in F, 



v 



Wz> + /CP (10) 



where K = V 3x6 ~ (0, F)V 6x6 and P = (N u N 2l N 3l N u N 2 , N 3 ). Then the charged- current 
Lagrangian in the mass basis is given by 

C cc = -^l L YuW- + h.c.^-^=l L ^(m + JCP)W-+h.c. (11) 

This mixing between the doublet and singlet components in the charged-current sector 
has several important phenomenological consequences, as listed below: 
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1. The flavor and mass eigenstates of the LH neutrinos are now connected by a non- 
unitary mixing matrix M = (1 — r))U, where the non-unitarity effects entering dif- 
ferent neutrino oscillation channels are measured by the parameter 77. In particular, 
the CP- violating effects in the leptonic sector will now be governed by the PMNS 
matrix Af instead of U through the Jarlskog invariant 



J% = Im (XJWV^) (12) 

where the indices a 7^ (3 run over e, \i and r, while % 7^ j can be 1, 2 and 3. In the 
standard PMNS parametrization of U by the three mixing angles 9y and the Dirac 
CP-phase 5, one can expand Eq. (fl2"|) up to second order in r) a p and Si 3 = sin6* 13 
(assuming those to be small) to obtain 

J&^J + Ajg,, (13) 

where the first term governs the CP- violating effects in the unitary limit and the 
second term gives the contribution coming from the non-unitarity effect: 

J = Ci2Ci 3 C23Si2Si 3 S23Sin5, (14) 



+ v^u^u;^ + v^u^u^u;,) (is) 

Note that the unitary term J vanishes if either s 13 — ► or 5 —>■ 0. However, 
AJ^ depends on the off-diagonal elements of rj (generally complex) and does not 
necessarily vanish even if both si 3 and 5 are zero; in fact, it might even dominate 
the CP- violating effects in the leptonic sector. 

2. The heavy neutrinos iVj and TVj entering the charged- current sector can also mediate 
the rare lepton decays, l~ — > l^j. Hence, unlike in the canonical seesaw model where 
this contribution is suppressed by the light neutrino masses [T7|, in this case it is 
constrained mainly by the ratio F = M^M^ 1 . The LFV decays mediated by these 
heavy neutrinos have branching ratios [10J 



BR(Z Q -> Im) 



a w s w m i a 
2567r 2 M^r 
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m 



(16) 



where T a is the total decay width of l a and the function I(x) is defined by 

2x 3 + 5x 2 — x 3x 3 In x 
[X) = 4(1 -x) 3 2(1 -xY { ' 

For degenerate RH neutrino masses, a reasonable assumption inspired by resonant 
leptogenesis [18], the amplitude is proportional to (fCK/ ) ~ \ FF^) , and hence, 

V / a/3 \ ) a(3 

for sizeable F and TeV-scale RH sector, one could expect appreciable rates in the 
LFV channels. On the other hand, in the conventional type I seesaw model, one has 
approximately /C/C^ = O^m^M^ 1 ^, and therefore, the branching ratio, BR(/ Q — > 
Ipj) oc O (ml) is strongly suppressed. 
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3. The heavy neutrinos iVj and iVj also couple to the gauge sector of the SM and can 
be produced on-shell, if kinematically accessible, at hadron colliders via the gauge 
boson exchange diagrams. Due to their pseudo-Dirac nature, the striking lepton 
number violating LHC signature of the fine-tuned type I and type III scenarios, 
namely pp — > ^/?+ jets, will be suppressed for heavy Majorana states due to 
cancellation between the graphs with internal lines of the N and N type which have 
opposite CP-quantum numbers. However, the LFV processes are insensitive to this 
effect and one can expect to get observable signals at the LHC. The most distinctive 
signature would be the observation of LFV processes involving three charged leptons 
in the final state plus missing energy, i.e. pp — > l^l^l^u(v)+ jets [19]. 

Thus we see that the phenomenology of the inverse seesaw mechanism depends 
crucially on the mass ratio F = MpM^ 1 . As noted earlier, we can choose the RH neutrino 
masses to be degenerate (with eigenvalue m N ), inspired by resonant leptogenesis. So we 
are left with a single mass parameter m^, together with the Dirac mass matrix and 
the arbitrary small mass parameter fi. In what follows, we explicitly determine the form of 
Md in the context of a realistic supersymmetric SO (10) model and then use the present 
experimental bounds on the elements of the non- unitary parameter \r]\ to get a lower 
bound on the RH neutrino mass scale mjv- Finally we fit the observed LH neutrino mass 
and mixing parameters by the inverse seesaw formula to determine the structure of fi. We 
then study the phenomenological consequences of our results. 



4 Embedding Inverse seesaw in realistic 50(10) GUT 

As we have mentioned earlier, in order to embed the inverse seesaw mechanism into a 
supersymmetric SO (10) theory, we have to break the B — L symmetry by using a 16 © 16 
pair rather than a 126 © 126 pair of Higgs representation. In this context, there are two 
symmetry breaking chains that are particularly interesting: 

. S o(10) ^ 3 C 2 L 2 R 1 B _ L ^ 3 c 2 L ly(MSSM) A ^ Y 3 c 2 L ly(SM) ^ 3c1q m 

• 50(10) ^ 3 c 2 l 2 r 1 B -l 3 c 2 L l /3fl l B _ L 3 c 2 L l y (MSSM) 

^ Y 3 c 2 L ly(SM)^3 c l Q 121] 

where, as an example of our notation, 3 C means SU(3) C . In this paper, we consider only 
the former (and simpler) case of 5*0(10) breaking chain. It was shown in Ref. [20J that 
it is possible to obtain the gauge coupling unification in this model with a low-energy 
(TeV-scale ) SU(2)r symmetry breaking scale Mr. However, they considered only one 
10h Higgs field which contains only a single bi-doublet (corresponding to the (1,2,2,0) 
representation of 3 C 2^2^1^_^). Getting a realistic fermion mass spectrum in this model 
is difficult (see, however, some recent ideas [22] on how this could be done). Instead, 
we consider a model with two 10 h at the TeV scale. This requires that we reexamine 
the unification issue with two Higgs bi-doublets. We show that we not only obtain the 
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gauge coupling unification at a scale consistent with the proton decay bounds, but also 
successfully reproduce the observed fermion masses and mixing. 

To study the running of the gauge couplings and the possibility of their unifi- 
cation at a scale Mq ~ 10 16 GeV, we divide the whole energy range (Mz, Mq) into three 
parts, according to the above mentioned symmetry breaking chain: 

• First, we have the well known SM from the weak scale Mz to the SUSY-breaking 
scale Msusy (which, for practical purposes, we assume to be a little higher than 

M z )\ 

• Then we have the MSSM from Msusy to the B — L breaking scale Mr (which is 
assumed to be of order TeV, so that it is of interest for colliders); 

• Finally, we have the Supersymmetric Left-Right (SUSYLR) model from Mr to the 
unification scale M G (expected to be around 10 16 GeV). 



equation 



The running of the gauge couplings at one-loop level is determined by the RG 



dai h 2 _ x _ _ x N ( bi (p,2 



_ —a { , or, a t {fi x ) = a, t (/x 2 ) + — In — 
am/2 Air Z7T \fj>i/ 



where aj = I 1 , p, is the energy scale, and b^s are the coefficients of the one-loop (3- 
functions. The SM and MSSM /3-functions are well known 



e = Q, ~ -7) , and bf SSM = (f , 1, -3) , (19) 

where i stands for ly, 2 L , and 3 C respectively. Before calculating the /3-f unctions for the 
SUSYLR model, let us first discuss the particle contents of this model. 

4.1 Particle Content of the SUSYLR Model 

Here we consider only the doublet implementation of the SUSYLR model [24], i.e. we 
use SU(2) doublets of the 16# Higgs field to break the B — L symmetry. In order to 
keep the model general, we allow for an arbitrary number of these doublet fields, to be 
denoted by til and hr respectively for SU(2)l and SU (2)r doublets. Likewise we have n w 
bi-doublets of the 10# Higgs field which, on acquiring VEVs, give masses to the fermions 
through Yukawa couplings. We also allow for an arbitrary number ns of singlet fields S a . 
This is the minimal set of particles in a generic SUSYLR model. 

However, it turns out that with this minimal set of particles, it is not possible 
to obtain the gauge coupling unification at a scale higher than ~ 10 15 GeV as required 
from current bounds on proton decay lifetime, r v ^ 10 34 years [25]. As we have shown 



'v 

later in Section 4.2, unification is possible only after adding the contribution from the 
color triplets ^3, 1, |j + c.c. which come from the 45# representation of Higgs at the 
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Multiplet 


SU(3) C 


SU(2) L 


SU(2) R 


U(1) B -l 


Q 


3 


2 


1 


+ 1/3 


Q c 


3 


1 


2 


-1/3 


L 


1 


2 


1 


-1 


L c 


1 


1 


2 


+1 


Xp 


1 


2 


1 


+1 


ryC 

Aq 


1 


1 


2 


-1 


Xp 


1 


2 


1 


-1 


v c 

A q 


1 


1 


2 


+ 1 




1 


2 


2 







1 


1 


1 





5 


3 


1 


1 


+4/3 


5 


3 


1 


1 


-4/3 



Table 1: The representations of the particles under the 3 c 2l2r1b-l gauge group in the 
doublet SUSYLR model. Here a = 1, nio, p = 1, ■■■,ni, q = 1, and a = 1, n^. 

The B — L quantum numbers given here are not GUT renormalized; to do so, we multiply 
a factor of */| (not y| as mentioned in Ref. [24]). 



unification scale Mg. It is justified in Appendix A that it is indeed possible to have these 
color triplets at TeV scale while all the other Higgs multiplets are still naturally heavy at 
the GUT-scale. 

The particle content and their representations under the Z c 2l2bXb~l gauge 
group are summarized in Table-1. Following the notation of Ref. [24], the SU{2) doublets 
and bi-doublets are represented as 

Other doublet pairs can be written in a similar way as the (Q, Q c ) pair. The charges of 
the fields must obey the relation 

Q = I 3L+ I 3R + ^—Ii (20) 
4.2 Gauge Coupling Unification 

The /5-function for a general supersymmetric model is given by [23] 

b s ^ SY = 2n g -3N + T(S N ) (21) 

for n g generations of fermions, the gauge group SU(N), and the complex Higgs represen- 
tation parametrized by T(Sn). For £7(1) gauge group, iV = in Eq. fT2TT) and the gauge 
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coupling is normalized as usual. For the particle content given by Table-1, the Higgs 
contributions in our SUSYLR model are explicitly given by 

3 

T 2L = nio + n L , T 2R = n w + n R , T 3c = 1, and T B _ L = A + -(n L + n R ) (22) 

Hence for three fermion generations, we find the /5-functions for our SUSYLR model to 
be 

6 susylr = + 3^ + + ^ n ^ + ^ _ ^ ^ (23) 

where i stands for 1 R -l, 2 £ , 2 R and 3 C respectively. Using these /3-functions, we can 
now obtain the running of gauge couplings up to the scale Mq, knowing the initial values 
at p, = M z [2B] 

a 1Y (M z ) = 0.016829 ±0.000017 

a 2L (M z ) = 0.033493^;~ 
a 3c (M z ) = 0.118 ±0.003 

and the matching condition [27] at jx = M R where the C/(l)y-gauge coupling gets merged 
wkoSU(2) R xU(l) B - L : 

a^(M R ) = \*&{M R ) + \a- B l _ L {M R ) (24) 

For numerical purposes, we assume Msusy = 300 GeV and M R = 1 TeV. 
Also we take the number of Higgs bi-doublets, nio = 2. However, the number of Higgs 
doublets can be arbitrary and we vary these parameters to get the unification. As shown 
in Figure- 1, we achieve the gauge unification for = and n R = 2, with the unification 
scale parameters 

Mq - 4x 10 16 GeV, and a^(M G ) ~ 20.3 (25) 

Note the asymmetry between ni and n R . We show in the appendix that since 
the VEV of the 45# Higgs breaks D-parity and decouples it from the SU(2) R breaking 
scale [28], it is possible to have only the right-handed doublets and no left-handed ones 
below the GUT scale. This leads to the asymmetry between a 2 L and a 2R , with ^£ ~ 1.3 
in our case. 



5 RG evolution of the fermion masses and mixing 

The RG evolution of the fermion masses and mixing have been extensively studied for 
both the SM and the MSSM cases [29J, but not for the SUSYLR model, even though the 
analytical expressions for the Yukawa couplings have already been derived in Ref. [24]. 
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Figure 1: Gauge coupling unification in the SUSYLR model. We have used n 10 = 2, n L = 
0, ur = 2, Msusy = 300 GeV and Mr = 1 TeV. As the running behavior is mostly 
controlled by the SUSYLR sector, the values of Msusy and Mr can be relaxed a little 
bit, still preserving unification. However, it should be emphasized that the choice of the 
number of bi-doublets and doublets is the only possible choice consistent with both gauge 
coupling unification as well as realistic fermion masses. Increasing uiq or n^, or changing 
riR in either way, will spoil the unification, and as already noted, reducing riio will not 
give us a realistic fermion mass spectrum. 

Here we present a detailed RG analysis in our SUSYLR model and obtain the quark and 
lepton masses and the CKM matrix elements at the unification scale Mq. 

The superpotential relevant for the RG evolution of the Yukawa couplings in 
the SUSYLR model is given by [21] 

W D ih a Q T T 2 $ a Q C + iti a L T T 2 $ a L C + iXapqXp^aXq + ^KpqXp^ 'aY \ 

+/M* ab S a Tt (Cr 2 $ 6 r 2 ) + i^ ap S a L T r 2Xp + i^ c q S a L cT r 2X c q (26) 

where we have suppressed the generational and SU(2) indices. Also we have ignored 
all non-renormalizable terms in the superpotential as their contributions to the RGEs 
are suppressed by Mr/Mq- We note that the superpotential given by Eq. (1261) has two 
additional terms of the form SL\ and SL C \ C (as required by the inverse seesaw model) 
as compared to that given in Ref. [21]. Also note that since the 5, 5 fields do not couple 
to any of the matter fields, they do not affect the renormalization group running except 
through their effect on the color gauge coupling evolution. 

We have seen from the previous section that the gauge coupling unification 
requires that we take two SU(2)r doublets and no SU{2)l doublet of Higgs fields. Hence, 
dropping the Xi X terms from the superpotential of Eq. ( |26l) . we have 

W D ^Q^^ + ^L^^ (27) 
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where a = 1,2; q — 1,2 and a = 1,2,3 corresponding to the two bi-doublets, two RH 
doublets and three fermion singlets, respectively. 

The renormalization group equations (RGEs) for the Yukawa couplings h a and 
h' a in Eq. (I2T1) are given by (with t = In fl) 



16tH 



,dh a 
~~dT 

dK 
dt 



h„ 



2h lh - y^s 



q9b-l 



+h b [Tr (3/ij> a + + 2/iJ/i a + 4 (/if ^ 



6« 



+K 



2h!lh' b 



o 2 o 2 d 2 

6 92L 6 92R 2 ^S-L 



Tr (3hlh a + h' b %) + 2h'tti a + 4 (>*V*) 6a + (/i LC )%S5 ba 



(2? 



(29) 



where the repeated indices are summed over and a, b = 1,2; g = 1,2; and a = 1,2,3 
corresponding to the two Higgs bi-doublets, two SU(2) R doublets and three fermion sin- 
glets respectively. Note that we have an additional contribution to the RGE of the lepton 
Yukawa coupling h' a as compared to those given in Ref. [21] which comes from the S\ C L C 
term in the superpotential. Note also the presence of the hb terms in the second line in 
both the Yukawa runnings even for a ^ b, which are characteristic of left-right models 
and are absent in the case of MSSM, arising from the Higgs self energy effects. 

The fermion masses arise through the Yukawa couplings h a and h' a when the 
two Higgs bi-doublets $i 2 acquire VEVs. In general, a linear combination of hi and h 2 
will give masses to the up-type quarks, and similarly different linear combinations for the 
other masses. The dynamics of the superpotential can be chosen in such a way that the 
bi-doublets acquire VEVs in the following simple manner: 





(30) 



and we identify the ratio v u /vd = tan (3 (MSSM) with Jv% + v% = 246 GeV. For numerical 
purposes, we use tan (3 (MSSM) = 10. To obtain the RGEs for the mass matrices, we choose 
the most frequently used renormalization scheme [29] where the Yukawa couplings and 
the Higgs VEVs run separately. The RGEs for the Higgs VEVs are obtained from the 
gauge and scalar self-energy contributions: 



16tt 



I6r 



. dv u 
dt 

>dvd 
dt 



|#2L + ~ Tr (3/4^2 + h$ti 2 



3 2 ^ 2 

2^2L + 2^2/? 



Tr (3/iI/ia + h'lhi 



4(a£V;) u 



(31) 
(32) 



Using Eqs. (1281 1291) for h a , h' a and Eqs. (13T1 1321) for v u , Vd, we have derived the RGEs for 
the physical fermion masses and the quark mixing in our SUSYLR model in Appendix B. 
Using the initial values for the mass and mixing parameters at fx = Mz 



m u {M z 



2.33^5 MeV, 



m c {M z 



677±|f MeV, 



m t (M z 



181 ± 13 GeV, 
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m d (M z ) = 4.69+°;™ MeV, m s {M z ) = 93At\H MeV, m b {M z ) = 3.00 ± 0.11 GeV, 
m e (M z ) = 0.48684727 ± 0.00000014 MeV, 

m^Mz) = 102.75138 ± 0.00033 MeV, m T {M z ) = 1.746691°;™ GeV, 

and with the quark-sector mixing parameters 9\ 2 = 13.04° ± 0.05°, 9i3 = 0.201° ± 
0.011°, 9 23 = 2.38° ± 0.06° and 5 13 = 1.20 ± 0.08, 

C12C13 S12C13 s 13 e~ lSl3 

V C km{M z ) = \ -s 12 c 23 - c 12 s 23 s 13 e lSl3 c 12 c 23 - s 12 s 23 s 13 e lSl3 s 23 c 13 

s 12 s 23 — Ci2C 2 3Si3e l<513 — C12S23 — Si2C23Si3e* <5l£i C23C13 

0.9742 0.2256 0.0013 - 0.0033i 

-0.2255 - O.OOOlz 0.9734 0.0415 

0.0081 - 0.0032s -0.0407 - 0.0007i 0.9991 

and the SM and MSSM Yukawa RGEs [29J we numerically solve the SUSYLR RGEs given 
in Appendix B to obtain the running quark and lepton masses and the CKM matrix 
elements at the unification scale M G : 

m u {M G ) = 0.0017 GeV, m c (M G ) = 0.1910 GeV, m t (M G ) = 77.8035 GeV; 
m d (M G ) = 0.0013 GeV, m s (M G ) = 0.0263 GeV, m b (M G ) = 1.7092 GeV; 
m e (M G ) = 0.0004 GeV,m M (M G ) = 0.0911 GeV, m T (M G ) = 1.7096 GeV 
/ 0.9793 0.2023 + 0.0018z 0.0005 - 0.0057z 



^ckm(M g ) 



-0.2023 + 0.0016i 0.9791 0.0240 | (33) 

V 0.0044 - 0.0056i -0.0236 - 0.0013z 0.9997 



We also have a mild running for tan/? with tan (5{M G ) = 7 from tan /3(Mr) = 10. 

Figure-2 shows the running of the quark and charged lepton masses up to the 
unification scale M G . Note that we are able to generate the fermion mass spectrum at 
the GUT scale with 

4.1,4.3,44 (34) 
m d 3 

Figure-3 shows the running of the CKM elements involving the third generation. Note that 
in addition to the significant running for the third generation CKM elements V u ^ c ^ td, ts, 
we have a relatively milder running for the other elements as well [cf. Eq. (1331 ], even in 
the third-generation dominance approximation. This is a characteristic of the Left-Right 
model, in contrast with the MSSM case where in the third generation dominance, the first 
and second generation elements do not run at the one-loop level. 



6 The Dirac Mass for Neutrinos in a Specific 5O(10) 
Model 

As discussed in Section 2, in order to implement the inverse seesaw mechanism, we have 
to use the class of SO(10) models in which the B — L subgroup is broken by a 16# © 16// 
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Figure 2: Running of fermion masses to the GUT-scale in our model for MgusY = 300 
GeV and Mr = 1 TeV. Note the b — t unification which is a generic feature of GUT 
models. 
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Figure 3: Running of the CKM mixing elements involving third generation to the GUT- 
scale in our model for Msusy — 300 GeV and Mr = 1 TeV. The running of other CKM 
elements, being small, is not shown here. 

pair. We also need at least two 10// and a 45// to have a realistic fermion spectrum. With 
this minimum set of Higgs multiplets {10^, 16//, 16#, 45//}, several 5*0(10) models 
have been constructed [30] . All these models require various dimension-5 operators to get 
right fermion masses: in principle, they are also present in our model. However, most 
of them e.g. ^-16j 16^16// 16//, are suppressed by the factor ~ 10~ 15 as the 16// 
Higgs acquires only TeV-scale VEV. The only other dimension-5 operator that can make 



14 



significant contribution to fermion masses is -^yl6jl6jl0//45//; we assume its effects to 
be small in our model and keep the dimension six operator 



■^16 i 16 J -10 H 45 H 45H 



(35) 



This operator is suppressed only by ~ 10 -4 as the 45// acquires a VEV at the 

scale Mq and plays an important role in the fermion mass fitting given below. 

The fermion mass splitting is obtained by the completely antisymmetric com- 
bination of the operator given by the expression (I35I) . i.e. in the notation of Ref. [31] 



(^ + \B[T i T j T k T l T m ]A ij A kl ^ m \^ + ) 



(36) 



with B — ]J[ T M and [...] denoting the completely antisymmetric combination. Here 

/x=odd 

$ and A denote the 10// and 45// fields respectively. When the following VEVs are 
non-zero: 

($9,10> 7^ 0, (A 12 , 34 ,56> ^ 0, (37) 

this antisymmetric combination acts as an effective 126// operator which gives the mass 
relation m e = -3m^ and m v = —3m u due to the VEVs (Aij), while m u and ma are split 
in the usual manner by the two 10// VEVs, (</»g,io)- To obtain a realistic fermion mass 
spectrum, we construct the following model using the Higgs multiplets {10//, 45//, 54//}. 
The 5*0(10) symmetry breaking to 3 c 2l2r1b-l is obtained by a combination of the 45// 
and 54//, with the following VEVs in an SU(5) basis: 

(45) oc diag(a, a, a, 0, 0), 

(54) oc diag(2a, 2a, 2a, 2a, 2a, 2a, — 3a, — 3a, — 3a, — 3a) (38) 
In this model, the fermion mass matrices at the GUT-scale have the following form: 



M u = h u + f } M d = h d + f } M e = h 



3/, M D = h u -3f 



(39) 



where the h Ut d matrices come from the usual Yukawa terms /iyl6j 16^10// (10#) and the 
/ matrix comes from the 45// contribution given by the expression (|35l) . where we have 
assumed the same coupling for both the 10// fields. The tilde denotes the normalized 
couplings with mass dimensions where the VEVs have been absorbed. We know the nine 
eigenvalues of the quark and charged lepton mass matrices at the scale Mq from our RG 
analysis in Section 5; however, we have 18 unknowns (for 3 hermitian matrices) to fit 
into Eq. fl39|) . Hence a unique fit is not possible; we just give here one sample fit that is 
consistent with all the masses and mixing at the GUT scale obtained from the RGEs. 

We work in a basis in which the charged lepton mass matrix is diagonal, i.e. 



( m e (M G ) 






«V(M G ) 






GeV 
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This immediately implies from Eq. (j2HD that 

/',/,, -U,. Vi^j (40) 

For simplicity, let us choose the /-matrix to be diagonal. Then Eq. (|4TJ|) implies that hd 
is also a diagonal matrix. We also have the following relations: 

hd,aa + face = TTl a , hd,/3/3 ~ = flip (41) 

where m a = (rrid, m s , m&) are the eigenvalues of Md and mp = (m e , m M , m T ) the 
eigenvalues of M e . These six equations (141j) now fix the hd and / matrices completely: 

\ (m d - m e ) 

\{m a - «v) 

\(m b - m T ) 

I 2.25 x 10" 4 \ 

-0.0162 GeV, 

-0.0001 / 

/ \ {3m d + m e ) 



./ 



V 



h n 



\ (3m s + m„) 







0.0011 
0.0425 
1.7093 







| (3m b + m T 



GeV 



(42) 



The h u matrix can now be determined by fitting to M u which, in this basis, is given by 

M u = VbKMMf^KM 

/ 0.0120 0.0384 - 0.0103i 0.038 - 0.4433i \ 

0.0384 + 0.0103i 0.2280 1.8623 + 0.0002z GeV (43) 

V 0.038 + 0.4433i 1.8623 - 0.0002i 77.7569 / 



Then from Eq. (139]) the h u matrix is given by 

/ 0.0118 0.0384 - 0.0103i 

h u = 0.0384 + 0.0103i 0.2442 

V 0.038 + 0.4433i 1.8623 - 0.0002i 

Hence the Dirac neutrino mass matrix is given by 

/ 0.0111 0.0384 - 0.0103i 

M D =\ 0.0384 + 0.0103i 0.2928 

I 0.038 + 0.4433i 1.8623-0.0002? 



0.038 - 0.4433i 
1.8623 + 0.0002z ] G<A" 
77.757 



0.038 - 0.4433i 
1.8623 + 0.0002* i G<A" 
77.7573 



(44) 



(45) 



It may be noted here that even though the specific form of the Dirac neutrino mass matrix 
may depend on the choice of the particular basis we have chosen, the individual values 
of the matrix elements are more or less fixed by the up-type quark mass values, due to 
the mass relation (13"9"|) . and hence, do not depend on the basis so much. Therefore, all 
the predictions of the model that follow from the form of Mu given by Eq. (145]) will be 
independent of the initial choice of our basis, upto a few 



/o. 
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7 Non-unitarity effects in the lepton mixing matrix 



In this section we obtain the non-unitarity parameter 77 using the structure of the Dirac 
neutrino mass matrix obtained in Eq. f|45|) and discuss the phenomenological consequences 
of our results. 

7.1 Bounds on \r]\ 

As discussed in Section 3, the non-unitarity parameter is given by 

7] ~ ]^FF ] with F = M D M N l (46) 

For simplicity, choosing Mn to be diagonal, and motivated by resonant leptogenesis, 
assuming degenerate eigenvalues for Mjv equal to m^, we have 

V * ^T M ° M d (47) 

With the form of Mo derived in the last section after extrapolation to the weak scale, we 
can readily calculate the elements of rj: 



1 GeV 



ni 



N 



2 f 0.1 0.0412 - 0.4144i 1.5134 - 17.247? 

0.0412 + 0.4144i 1.78 72.6794 - 0.0005i j (4? 

\ 1.5134 + 17.247i 72.6794 + 0.0005i 3024.93 



This is to be compared with the present bounds on \rjij\ (at the 90% C.L.) 





( 2.0 X 10~ 3 


3.5 x 10~ 5 


8.0 x 10- 


-3 


M < 


3.5 x 10~ 5 


8.0 x 10~ 4 


5.1 x 1(T 


-3 




V 8.0 x 10~ 3 


5.1 x 10~ 3 


2.7 x 10- 


-3 



(49) 

10~ 3 J 

This gives a lower bound on the mass of the RH neutrino: 

m N > 1.06 TeV, (50) 

which should be kinematically accessible at the LHC to be produced on-shell. Note that 
the right handed neutrinos are pseudo-Dirac fermions in our model (with small Majorana 
component) which is distinct from the type I seesaw models where they are pure Majorana. 
As a result the like sign dilepton final states which are the "smoking gun" collider signals 
of type I seesaw are suppressed in our model; however, the trilepton signals can be used 
in this case for testing these models [T9~] . 

With this lower bound on m^, we get the following improved bounds on \ijap\- 

\r) ee \ < 8.9 x HT 8 , \rj efl \ < 3.7 x 10" 7 , |r/ eT | < 1.5 x 10~ 5 , 

\rjw\ < 1.6 x 10~ 6 , \ VlMT \ < 6.5 x 10" 5 (51) 
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At least one of these bounds, namely \r) efi \, is reachable at future neutrino factories from 
the improved branching ratio of fi — > ej down to 10~ 18 [33]. Similar sensitivities are 
also reachable in the PRISM/PRIME project [33]. We note that relaxing the condition 
of degenerate RH neutrinos but fitting the neutrino masses affects the values of r\ a p ; we 
present these results in Table-2. It appears that \r] efl \ values are all accessible to the 
future /i — > e + 7 searches; The largest value of \rj^ T \ in this table may also be accessible 
to neutrino oscillation experiments, preferably with short baseline (L ^ 100 km). 



itin 2 mN 3 


M \Ver\ KA 


1100 1100 1100 
100 100 1100 
50 50 1200 
30 30 2100 


3.7 x 10~ 7 1.5 x 10~ 5 6.5 x 10~ 5 
7.9 x 10~ 7 1.6 x HT 5 8.9 x 10~ 5 
2.5 x 10~ 6 2.2 x HT 5 1.6 x 10" 4 
6.7 x 10" 6 4.4 x 10" 5 3.2 x 10" 4 



Table 2: predictionss for the non-unitarity parameter |r] Q/ g| for the above choice of param- 
eters in the model including RH neutrino masses (given in GeVs). 

7.2 Fitting the Neutrino Oscillation Data 

The structure of the small mass parameter /1 can be obtained using the inverse seesaw 
formula, Eq. ([1]): 

/1 = F~ x m v [f tS ) X (52) 

where m v is diagonalized by the new PMNS matrix M = (l — r])U instead of U in Eq. ([6]): 

m v = Mm u M T (53) 

The form of U is obtained from the standard PMNS parametrization using the 2a results 
from neutrino oscillation data [35J 

Ami = 7.67 (l^fr) xlO" 5 eV 2 , 

Am 2 tm = 2.39 (l™) x 10- 3 eV 2 , 

sin 2 £ 12 = 0.312 (lj^jg) , 

sin 2 £ 23 = 0.466 (l^; 2 ? 2 ) , (54) 

Here we assume #13 = 0. Now using the form of rj obtained in Eq. fj4*8l) and taking 
rriN = 1.1 TeV for its lower bound value, we get the new PMNS matrix M = (1 — f])U. 
For illustration, let us assume normal hierarchy for neutrino masses with mi = 10~ 3 eV. 
Then we obtain from Eq. ( |52l ) 

/ -1.5934 + 0.0283* 0.2244 - 0.0063i -0.0044 + 0.0092i \ 
y. ~ 0.2244 - 0.0063i -0.0322 + 0.0012i 0.0006 - 0.0013i GeV (55) 

\ -0.0044 + 0.0092i 0.0006 - 0.0013i 4.0 x 10" 5 + 5.1 x 10 5 i / 
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7.3 CP-violation effects 



The CP-violation effects due to non-unitarity are measured by the Jarlskog invariant 
AJ^j given by Eq. (fl"3l) . Note that AVg is non-zero in our case as r\ is a complex matrix 
(the phases arising from the Dirac neutrino sector). Using the values of obtained from 
neutrino oscillation data given by Eqs. (1541) and the structure of r] determined in Eq. (I48p 
with m^v = 1-1 TeV, we obtain the following values for AJ^: 

AJ e 12 ~ -2.4 x 10~ 6 , (56) 

AJ% ~ -2.7 xlO" 6 , (57) 

Ajf T ~ 2.7 x 1(T 6 , (58) 

Ajf T ~ 2.7 x 10~ 6 , (59) 

AJ r 12 ~ 7.1 x 10~ 6 (60) 

and AJ^ = AJ^ = — AJ^ 2 = AJ 2 f = AJ^. Note that these values are just one order 

of magnitude smaller than the quark sector value, Jqkm = (3.05ioio) x 1CT 5 [26] . and 
can be the dominant source of CP- violation in the leptonic sector for vanishing # 13 , thus 
leading to distinctive CP- violating effects in neutrino oscillations [361 [37]. F° r instance, 
the transition probability for the "golden channel" z/ M — > v T with non-unitarity effects is 
given by 



l9 9 9 9 /Amo 1 L\ „ ( Am\,L\ , s 

P^r ^ Mv^\ + 4s 2 3 c 2 3 sin 2 I ^ I -4|^ T |sin^ T g 2 3C 23 sin I ^ I (61) 

where the last term is CP-odd due to the phase <5 MT of the element rj^ which, in our 
model, is ~ 7 x 10~ 6 [cf. Eq. (14"81)]. Hence, the CP- violating effects should be pronounced 
for long-baseline neutrino factories. 



7.4 LFV decay rates 

Lepton flavor violating decays such as /i — > ej, r — > and r — > fij are often a signature 
of seesaw models for neutrino masses. In this model, they can arise from the non-unitarity 
effects and can be obtained using Eq. ([To]) which, for degenerate RH neutrinos, becomes : 

2 



BR(l a Ipy) 



a W s W m l a 

2567r 2 M^r c 



(iCK)) I 

V J a0 



in 



N 



(62) 



with /C = V3X6 and I(x) defined in Eq. ( Tl7l) . Now that we know all the three 3x3 mass 
matrices entering the inverse seesaw formula given by Eq. ([3]), we can easily determine 
the structure of the full unitary matrix V by diagonalizing the 9x9 neutrino mass matrix 
Ait,, and hence, obtain 

The total decay width T a entering Eq. (1621) is given by h/r a where the mean 
life for /i and r are, respectively [26] 

= (2.197019 ± 0.000021) x 10" 6 sec, 
r T = (290.6 ± 1.0) x 10~ 15 sec. 
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Using these values, we obtain the following branching ratios for the rare LFV decays 



BR(/x - 


-> ej) 


~ 3.5 x 10~ 16 , 


(63) 


BR(r - 


-> ery) 


~ 1.1 x 10~ 13 , 


(64) 


BR(r - 


+ Hi) 


~ 2.0 x 10~ 12 


(65) 



We have estimated the contribution to /x — > e + 7 branching ratio from the off diagonal 
Dirac Yukawa coupling contribution to slepton masses and find that for universal scalar 
mass of 500 GeV and tan (3 ~ 5, it is comparable to this value or less. Such values for 
/x — > e7 branching ratio are accessible to future experiments [33, 34J capable of reaching 
sensitivities down to 10~ 18 . They can be used to test the model. 

In our model we assume that squark and slepton masses are above a TeV so 
that their contribution to the flavor changing neutral current effects are negligible. The 
predictions for \x — > 3e and /x — > e conversion [38] for a TeV-scale slepton mass, as in our 
model, are much smaller than what can be probed in planned experiments. 

8 Summary 

In conclusion, we have presented a TeV scale realistic inverse seesaw scenario that arises 
from a supersymmetric 5*0(10) model consistent with gauge coupling unification and 
fermion mass spectrum. This required us to carry out an extrapolation of quark masses 
and mixing to the GUT scale with a TeV scale SUSYLR rather than MSSM. This appears 
to be the first time that such an extrapolation is carried out. Implementation of inverse 
seesaw within the 50(10) helps to reduce the number of parameters making the model 
predictive. We present our expectations for the non-unitarity of the PMNS leptonic mix- 
ing matrix with the choice of parameters and its other phenomenological consequences. 
The heavy RH neutrinos which are pseudo-Dirac fermions have TeV scale mass and can be 
produced in colliders, thus giving rise to distinctive signatures. We also give our predic- 
tions (with our choice of parameters) for the non-unitarity contribution to the branching 
ratios for the rare LFV decays of muons and taus. The model can also be tested by the 
production of Wr and Z' bosons which are at the TeV scale. Of these, the branching 
ratio /x — > e + 7 could be testable in future experiments. Some of the elements of the 
non-unitarity matrix \rj\ predicted by our model may be accessible to the next generation 
neutrino factories too. 
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Appendix A: Masses of the 50(10) Higgs multiplets 



As discussed in Section 4.2, we obtain the gauge coupling unification at an acceptable scale 
only after including the contribution from the color triplets 5 (3, 1, 1, +|) , 5 (3, 1, 1, — |J . 
This pair of Higgs fields is contained in the 45 representation of Higgs in a generic SO (10) 
model. However, in principle, there could be other light gauge multiplets of 45 and/or 
54 that might contribute to the gauge coupling running as well. Here we argue that in 
a generic SO (10) model with only 45 h and 54# representations of Higgs (apart from 
the essential 10 H and 16#), it is possible to have only the S's as light states (TeV scale) 
whereas all the other states are very heavy at GUT scale, and hence, do not contribute 
to the RG running. It turns out that we need to have at least two 45#'s in our model in 
order to have these light color triplets. 

The most general Higgs superpotential with two A = 45's and a E = 54 Higgs 
fields is given by 

W H = ^miA 2 + ^A' 2 + ^m 2 E 2 + AiE 3 + A 2 EA 2 + A 2 EA' 2 + A 3 EAA' (66) 

where we have absorbed the AA' term by a redefinition of the fields. The Higgs fields 
A, A' and E contain three directions of singlets (with A and A' VEVs parallel) under 
the SM subgroup 3 c 2l1y [39J. The corresponding VEVs are defined by 

2 2 

(A) = J2AA, (A'HE^'A (E) = EE (67) 

i=l i=l 

where in the notation of Ref. [39], the unit directions A\ and E in the Y-diagonal basis 
are given by 

90] = A' i, 
+ 34 + 56] = A' 2 , 

-2 x [12 + 34 + 56] + 3 x [78 + 90]) (68) 

where the upper and lower indices denote the 3 c 2ily and 4 c 2l2r quantum numbers 
respectively. The unit directions in Eq. (|67|) satisfy the orthonormality relations 

A\ ■ Aj = 5ij and E ■ E = 1 (69) 

The superpotential of Eq. (166]) calculated at the VEVs in Eq. fl67|) is given by 



M 


4(1.1,0) 
— ^(1,1,3) — 


1 r 

,[78 


A 2 


4(1.L0) 
- A (15,l,l) " 




E 


£5(1,1,0) 
~ ^(1,1,1) ~ 


1 

Veo 



(Wh) = l -m x (Af + ~m' 1 (AY +\m 2 (E)* 

+\ l (Ef + X 2 (E)(A) 2 + X' 2 (E)(A') 2 + X 3 (E)(A)(A') 
= \mr(A\ + A\) + ^[(A? + A'i) + \m 2 E* + ^=E* 



L \ 2 (3Aj - 2A 2 2 ) + \' 2 (3A'? - 2A' 2 2 ) + X^A^ - 2A 2 A' 2 ) (70) 



2V15 
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using the definitions in Eqs. (|68|) and the orthonormality relations given by Eqs. 
The VEVs are determined by the minimization of the superpotential with respect to the 
fields: 







d d d d 



(W H ) = 



[dAi dA 2 ' Mi' dA 2 ' dE 
This yields a set of five equations for Ai, A 2 , A[, A 2 and E: 

= m x A x - 



(71) 



3 X 2 EA 1 + -4=A 3 £Ai, 










m 1 A 2 - 
m' 1 A 1 + 
m[A' 2 - 
m 2 E + 



15 
2 



15 



XoEAo 



15 
2 



\' 2 EA l + 



15 



XnEA' 



2 






2 


2 


s/l5 




3 


2 


s/15 




2 


2 





2" 



X 3 EA 1 , 
X 3 EA 2 , 



(72) 



2V15 



3Ai£ 2 + A 2 (3A? - 2^) + X' 2 (3Af - 2A 2 2 ) + X 3 (3A 1 A' 1 - 2A 2 A' 2 



As in our model, the SO (10) symmetry is broken by the 45 and 54 VEVs to Z c 2l2bXb-l 
gauge group at the scale Mq, we are interested in the 3 c 2l2r1b~l symmetry solutions 

A, = A[ = 0, A 2 + 0, A' 2 + 0, E^O 

Hence it follows from Eqs. (1721) that 



2X 2 E X 3 E A' 



mi 



15 



15 A 2 ' 



nil 



2X' 2 E _ X 3 E A 2 

~JEa f 2 



15 



(73) 



In order to study the mass matrices, it is convenient to decompose the Higgs 
representations under the SM gauge group 3 C 2 L 1 Y - In Table-3 we present the explicit 
decompositions of all the Higgs representations under the chain of subgroups 

4 c 2^2 j r D 3 c 2l2r1b-l D 3 c 2jAy- 

Using the Clebsch-Gordan coefficients given in Ref. [39], we obtain the masses of these 
multiplets as follows. The basis designating the columns (c) of the mass matrices is given 
in the same way as in Table-3 while the rows (r) are designated by the corresponding 
complex conjugated 3 c 2i\y multiplets. 



First, we obtain the masses of the multiplet ^3, 1, |) + c.c. 



in the basis 



c . ^(3,1,1). . A , 

L • ^(15,1,1)' A (15,1,1)' 1 ' ^(15,1,1) ' 



(3,1,-|) 



mi — 



2X 2 E 



15 



X3E 




(15,1,1) 

XzE 



-A 



-1 

M 
A> 



(74) 
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using Eq. (1731) . It is obvious that det(Ms) 
zero while the other eigenvalue is given by 



0, and hence, one of the two eigenvalues is 



Tr(M s 




(75) 



The zero eigenvalues (six in total) are easily identified as the longitudinal Nambu-Goldstone 
modes as the SU(4) C gauge group breaks to SU(3) C x U(1)b-l an d they acquire mass of 
order M G by the usual Higgs mechanism once the 45# gets VEV at the GUT scale. We 
keep the other six eigenvalues given by Eq. (!75|) at TeV scale by fine-tuning the coupling 
A3. In what follows, we explicitly calculate the mass eigenvalues for all the other multi- 
plets given by Table-3 and show that it is possible to have only the above six massive 5's 
at the TeV scale while all the other states of 45 and 54 are heavy at the GUT-scale. 

We note that once we assume A3 to be small, the effect of the second 45// mul- 
tiplet becomes negligible and we can as well drop the primed terms in the superpotential. 
For simplicity, we also assume that A 2 = E ~ Mq. Then the VEV conditions given by 
Eqs. (J72D yield 



777.1 



2\ 2 E 



E 



m 2 



A, 



. 4(1,1,0) 2(1,1,0) 

■ (1 1 q-1 1 



p(l,l,0) 

(15,1,1)' ^(1,1,1) > 1 ' -^(1,1,3) J 
2X 2 E 



l (l,l,3)> 



(1,1,0) 2(1,1,0) ^(1,1,0) 
^(15,1,1) J "^(1,1,1) 



mi — 



V15 
2\ 2 A 2 



\ 



3A 2 Ai 

VTE 



2X 2 A 2 



15 



15 



'lb 

1 1 1 3Xi E 

m2 + w / 







-2A 2 



So the mass eigenvalues are 



Mj 



(1,1,0) 



5E 



M. 



(1,1,0) 



15 
E 



M, 



(1,1,0) 



2V15 
E 



2V15 



A 2 + 2 A i 



A 2 + -Ai 




A 2 + ^Ai) +16A1 



3 2 
2 



A 2 + ^A!) + 16A1 



(76) 



We list below the mass eigenvalues for all the multiplets given in Table-3. 
• (1,1,0) : We have three such states and the mass matrix is given by 





-2A 2 
A 2 + |Ai 



7^0, 

+ (77) 



[(1, 1, 2) + c.c.] : There is only one such multiplet and its mass is 



c ■ /4 (1 ' 1,2) r 
c • ^(l.i^)' I 

M^ 1 ' 2 ) = mi + 



X(l,l,-2) 
^(1,1,3) 

\iE~- 



5E 



15 



A 2 ^ 



(7f 
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SO(10) 


4c,2 L ,2 fl 


3 C , 2j,, 1_b-l 


3 C , 2 L , ly 




(1,2,2) 


(1,2,2,0) 


(1,2,±1) 


10 


(6,1,1) 


(3,1,1,-1) 
(3,1,1,1) 


(3, 1,-1) 
(3,1,1) 




(4,2,1) 


(3,2,1,|) 
(1,2, 1, -1) 


(3,2,|) 
(1,2, -1) 


1U 


(4,1,2) 




(3 1 1) 
[0, 1, 3 j 

(3,1,-!) 






(1,1,2,1) 


(1,1,2) 
(1,1,0 




(1,1,3) 


(1,1,3,0) 


(1,1,2) 

(1,1,0) 

(1,1,-2) 




(1,3,1) 


(1,3,1,0) 


(1,3,0) 


45 


(6,2,2) 


(3,2,2,-1) 


(3,2,|) 
(3,2,-|) 






(3,2,2,1) 


(3,2,|) 
(3,2,-1) 




(15,1,1) 


(1,1,1,0) 

(3,1,1,|) 
(3, 1,1,-|) 
(8,1,1,0) 


(1,1,0) 

(3,14) 

(3,1,-1) 
(8,1,0) 




(1,1,1) 


(1,1,1,0) 


(1,1,0) 




(1,3,3) 


(1,3,3,0) 


(1,3,2) 
(1,3,0) 
(1,3,-2) 


54 


(6,2,2) 


fqo 9 _2\ 
V 3 '^' 3 J 


(<3,A 3 J 

(3,2,-!) 






(3,2,2,1) 


(3,2,|) 
(3,2,-|) 






(6,1,1,-1) 


(6,1,-!) 




(20', 1,1) 


(6,1,1,!) 

(8,1,1,0) 


(6,1,!) 

(8,1,0) 



Table 3: Decomposition of the 10, 16, 45 and 54 Higgs representations under the chain 
of subgroups 4 C 2 L 2 R D $ C 2 L 2 R 1 B _ L D 3 c 2 L l y . 
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(3, 2, — |) + c.c. : There are two such multiplets and the mass matrix is 

3-2,-f). 



L ■ ^(6,2,2) 



E) 



mi ' 2V15 
A 2 Ai _ A 2 A 2 

2 VE 



(6,2,2) 
A 2 _E A 2 Ai 



. _ . 4(3.2,1) ~(3,2,§) 
' 1 • (6,2,2) ' (6,2,2) 



A 2 A 2 

2 V6 
3Ai E 



m 2 + 



2a/15 



with the eigenvalues 



EX 2 
2VTE 



4± V14 



X2E 

VT5 \ -Jl 



7^0 



(3,2,|)+ c.c. 



: There are two of them and the mass matrix is 



(3,2,i N 
l (6,2,2) 



' ^(6,2,2) 



^(3,2,-1) -(3,2,-1) 
^(6,2,2) > ^(6,2,2) 



1 A 2 i? A 2 Ai 
'"1 o./Ts 2 



A 2 A 2 



2V15 

A 2 Ai A 2 A 2 1 3Aii? 

' 7T m 2 + 



XoE 



3 



2 . 7T ' ■ ' - ■> t : ; / v 1 5 

with the same eigenvalues as the previous one: 

EX2 



2V15 



4± V14 



7^0 



(1,3,0) : There are also two of them and the mass matrix is 



„ . 4(1.3,0) ft(l,3,0). . 4(1,3,0) p 
L • ^(W)' -^(1,3,3)' 1 • ^1.3.11' rj - 



m 1 + 



3A 2 g 



15 



A 2 ^i 



l (l,3,l)' 

E 



(1,3,0) 
(1,3,3) 

5A 2 



A 2 ^i 



m 2 + 



9XiE 

VTE 



So the mass eigenvalues are 



Mi 



(1,3,0) 



5£ 



A2 ^ 0, Mp<°) = -|= (A 2 + H Al ) + 



[(1, 3, 2) + c.c] : There is only one such multiplet whose eigenvalue is given 



6(1,3,2) 
(1,3,3)' 1 



E 



(1,3,-2) 
(1,3,3) 



M(1,3 ' 2) = ™2 + A ^ = "7= f A 2 + ^Ax) ^ 

(6, 1, — |) + c.c. : Its eigenvalue is 



( 6.1.-1) 



C : -^(20',l,l/> r ' -^(20', 1,1) 



m( 6>1 -5) 



m 2 - 



« A,E = +L (a 2 - £0 * 



unless A 2 = ^Ai (which we assume not to be the case). 
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(8,1,0) : There are two of them and the mass matrix is 



c : A 



(8,1,0) 



E 



(8,1,0) 



(15,1,1)' ^(20', 1,1)' 

mi 



r:A 



2X2M /2 XA 



\ 2 A 2 



with the mass eigenvalues 



Mi 



,1,0) 



E 



m 2 



(8,1,0) 
(15,1,1)' 

E 



E, 



± 



,1,0) 
(20', 1,1) 





15 \ ,/f A 2 



15 

y 



A 2 -^Ai) + 90A 2 , 



A, 



/45 -A 2 



7^0 



(84) 



Thus we see that all the other multiplets have non-zero masses, and moreover, all these 
masses are of order E ~ M G . Hence, none of these multiplets will contribute to the 
running of gauge coupling up to the unification scale Mq except the color triplets since 
these color triplets have masses of order of the SUSY breaking scale. 



3,1, 



c.c. 



Note that the 10^-Higgs field also has a color triplet pair 
under the SM gauge group, apart from the TeV-scale bi-doublet fields $i )2 used in the 
SUSYLR model in Section 4 which reduce to (1, 2, ±1) under the SM gauge group. At 
the GUT-scale, the H = 10# field interacts with the E = 54# field by the following term 
in the superpotential: 



W 



= im 3 H 2 + A 3 EH 2 



(85) 



After the 54// acquires a VEV, this gives rise to the color triplet mass 



^ (6,1,1) 



m( 3,1 ^) 



m 3 - 



2X 3 E 



(6,1,1) 



15 



(86) 



while the doublet mass is 



r . rr(l,2,l). . rfO-A-l) 
C • 11 (1,2,2) i r • n < 



'(1,2,2) 



M^ = m 3 + J^\ 3 E 



(87) 



We see that the (1,2, ±1) field can be made light by fine-tuning m 3 + J^X S E ~ TeV 



which still leaves the \3, 1, 3 J field heavy (of order Mq). 

Finally, let us discuss how only the right handed doublets fields (x c , X c ) from 
16//-Higgs fields (iPh) remain massless at the GUT scale. Note that in the left-right 
language, the fields in 16 are Q H (3,2, 1, §) © Q C H (3, 1,2, —J) and x (1,2,1,-1) © 
X°(l, 1, 2, +1), and similarly for 16 H = ip H field. The superpotential involving these 
fields is 

W w = MieTpHipH + \tp H Aip H 
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The second coupling has been worked out explicitly in Ref. [ID]. On substituting the 
VEV of the 45-Higgs field (A), we get the following masses for the Qh (3, 2, 1,|) © 
Q% (3, 1, 2, -§) and X (1, 2, 1, -1) © x c (l, 1, 2, +1) fields: 

M Q H -Q H = M i6 + XA 2; M Q* H -Q% = M i6 - XA 2 ; 

M x _^ = Mi 6 - 3XA 2 ; M x c_ r = M w + 3XA 2 (89) 

From this we see that to get only the \ c fields light, we have to fine-tune M%e + 3XA 2 ~ 
TeV. With this assumption, all other fields remain heavy at the GUT scale. 



Appendix B: RGEs for fermion masses and mixing 



Given the form of the bi-doublets VEVs as in Eq. (130]) . it immediately follows from the 
first two terms of the superpotential Eq. (1271) that the fermion mass matrices can be 
written as 

M u = -^u u h 2 , M d = -j=v d hx, M e = -^VdK, and M D = ^v u ti 2 (90) 

Henceforth, for clarity, we will denote the Yukawa couplings as 

hjj = h 2 , Yid = h\, h E = hi, hj\j = h 2 

Then using Eqs. ( 1281 [29]) and fl3T| [321) the RGEs for the fermion mass matrices can be 
written as 

,dM„ 



16?r 2 



1Qti z 



dt 



dt 



Ah) u h u + 2h ] D h D -Y J C\ q) g 2 



+M d tan(3 \Tt (3h^ D h v + h E h N ) + 2h ] D h v + C 



M, 



Atf D h D + 2tf u h u -Y^C\ q) gt 



tan (3 



Tr (3h\jh D + h ] N h E ) + 2h\ T h D + G 



16tt" 



,dM e 
~dT 



Ah E h E + 2h ] N h N + c*-Y. cf ] g 2 



M. 



D 



tan p 



Tr (3hth D + h ] N h E ) + 2hlh E + G 



■21 



16n 2 



,dM, 



D 



dt 



M, 



D 



Atf N h N + 2h E h E + - Cfg\ 



+M e tan/3 [Tr (3h ] D h v + h E h N ) + 2h E h N + Cf 2 



where C% = 4 (/i*V*) 



a 



t 

at; 

16 3 3 1 
~3~' 2'2' 6 



and for z = 3 C , 2 L , 2 K , 



> °i 
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} 3 3 3 

'2' 2' 2 



(91) 



(92) 



(93) 



(94) 



(95) 



Note that the second line in each of the above mass RGEs, Eqs. (19 1H94I) . is characteristic 
of the left-right models, and does not appear in MSSM. 

Not all the parameters of the Yukawa matrices are physical. Under an arbitrary 
unitary transformation on the left (right) -handed fermion fields, Fl(r) ~ > 
(where T = U,D, E, N), the Yukawa matrices undergo a bi-unitary transformation, hf — > 
LfhfR^ and the charged current becomes off-diagonal, with the CKM mixing matrix 

Lxjl} D . We will also have a leptonic counterpart of the CKM matrix that represents the 
mixing between the charged lepton and Dirac neutrino sector. However, as the running of 
lepton masses is very mild and we are working only to the one-loop order, we can safely 
ignore this mixing in the leptonic sector. Moreover, if we assume the CP phase in the 
Higgs VEV to be zero, then the mass matrices are Hermitian and Lf = Rf (manifest 
left-right). Thus we may perform scale-dependent unitary transformations £/(/i) on the 
fermion bases so as to diagonalize the Yukawa matrices, and hence the mass matrices, at 
each scale: 

hf{fx) = L f (fx)h f (fx)L\(fi), and M f = L f (fx)M f (fx)L\^), (96) 

where hf and Mf denote the diagonalized Yukawa and mass matrices, respectively. 

The RGEs for the physically relevant quantities, namely the mass eigenvalues 
Mf([i) and the scale-dependent CKM matrix Vckm(aO = Lu{^)L} D {^), are both contained 



in the RGEs of M](/i) = L}(ii)M f (ii)M}(fj,)L f (ii): 

1 



l v l\j,mZ 



16tt 2 



Ah 2 v + 2h 2 D 



{q) i 



2Mi 



167T 



tan/? [{Tr {W CKM h D V^ K Jiu) + C*} (V CKM M d V^ KM M u 



+2V C KMM d h D V ( l KM h u M u + h.c. 



(97) 



dt 



LdL^ d M 



+ 



16vr 2 
{Tr (3V C 



2Mi 



dt 
d_ 

dt 



1 1 

167r 2 tan j3 

+2M d h D V ( t KM h u M u VcKM + h.c. 
1 



+ C? 2 }(M d v£ KU M u V c 



CKM 



LeLe, Mi 



(M 



LnL^ n , M]\ 



+ 



+ 



167T 2 
1 

16^2 



4/i| + 2h\ + m. (C*) 

Ah 2 N 



i 

+ 2h 2 E + M(C*)-j:cPgl 



2Mi 



2M 2 D 



(98) 
(99) 

(100) 



where L 



§ and 3« (C7*) 



denotes the real part of C x . The commutator 



LfLf, Mf 



has vanishing diagonal elements because Mj is diagonal. Thus the RGEs for the mass 
eigenvalues mj follow immediately from the diagonal entries of Eqs. f)97H100p . Using 
dominance of Yukawa couplings of the third generation over the first two, i.e. 



hf > h 2 c > h 2 u , 



hi > hi > h 2 



d- 



» K » h 2 e , 



h% 3 »/4 ^>h 2 Nl , 
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we obtain the following RGEs for the mass eigenvalues of the fermions: 

2 d/TYliy^ 



167^ 
167T 2 "'"' 



dt 
■,dm^ 



16tt 
16tt ; 

167T 
167T 
167T 
167T 

16tt 
IQtt 
16tt 

167T 



,dmt 
~dT 

drrid 
~dT 

2 dm s 



dt 

2 dm b 
~dT 

2 dm e 
~df 
.dm, 



2 Lb " b H 

dt 
2 dm T 
~dT 

2 dm^ 1 

dt 
2 dmj\[ 2 
dt 

2 dfTljsf^ 

dt 



rn, 



m d + 



m, + 



m b + 



\ i 

Uh 2 c + 2h 2 s -^ ) g 2 

\ i 

U? + 2/^-£C^ 2 

\ i 
\ i 
\ i 

Lhl + 2hl-Y.C\ q) g> 

\ i 

Uhl + 2h% 1 +r l -Y,cPg> 

\ i 

Uhl + 2hl 2+ r l -Y.Cf ) g} 

\ i 

Uhl + 2h% 3 +r l -^CPg^ 

\ i 

Ul i+ 2hl + r l -j:c! ) g* 

\ i 

Uhl 2 +2hl + r l -Y,Cf ) gj 

\ i 

Uhl, i+ 2hl + r l -Y.Cf ) g* 

\ i 



m u + tan/3 3\V tb \ 2 h b h t + r q j ^ \V t 

j=d,s,b 



+ tznf3[3\V tb \ 2 h b h t + r q } ]T IK 

j=d,s,b 



cj\ m j 



m t + tan/3 [{3\V tb \ 2 + 2) h b h t + r q ] \V t 
1 



tb\ 2 m b 



tan (3 
1 

tan p 
1 

tan (3 
m e 

«V 
m T 

m Nl 

m N2 
m N „ 



3\V t b\ 2 h b h t + r q ] £ \ V id\ 



j=u,c,t 



3\V tb \ 2 h b h t + r q ] \ V J 



v js \ rrij 



j=u,c,t 



(3\V tb \ 2 + 2) h b h t + r q ] \V t 



tb\ 2 m t 



(101) 



where r q = 3? (Cf 2 ) and r, = 3? (C x ). 



The VEV RGEs, Eqs. (1321) and (l3Tj) . for third generation dominance become 

(102) 



, a 2 dv u 

l07T — — ~ V u 

dt 

l07T —7— ~ Vd 



dt 



2 _L 2 'i/i 2 /i 2 Z^* 

~2p2L < l^iiR t N 3 22 



(103) 



The RGE for the CKM matrix Vckm = LuL^ D is given by 

-jjVcKM = LjjLd + LjjD d = LuL\jVckm — VckmLdL\), 
or, —V a/ 3 = E (^^^ v tP~ V ai (L D L ] D ^ 



dt 



~f=U,C,t 



y=d,s,b 



7/3 



(104) 
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However, the diagonal elements of L UD L\j D are not determined by Eqs. (J97|) and 
This is because Eq. ( 195]) determines Z^d only up to right multiplication by a diagonal 
matrix of scale-dependent phases. These undetermined phases contribute arbitrary imag- 
inary functions to the diagonal elements of Lu,dL\j d . But the off-diagonal elements are 
unambiguously determined because they receive no contribution from the phases. We 
can, nevertheless, make the diagonal entries of L^D-^jy D , which are manifestly imaginary, 
vanish by an appropriate choice of phases. With this choice of phases, we can then obtain 
the RGEs for the CKM matrix elements using Eq. (11041) : 



d v 



r y=u,c,t 
7^q 



D 



E {^ uL u) ai v ip - E v ^(l d l[ 
l 



■y=d,s,b 



7/3 



16vr 2 



E 



. J=U,C,t 



tan P 



m- 



+ 



4 m 2 a + m 2 
v d m a~ m y V 



{Tr (Wh D V^h 
( VM 2 d V ] ) \ : 



u 



,} (VM d V^) 



°7 



/ 07 



7/3 



~ E Krr 



■y=d,s,b 
7^/3 



+ 



tan/3(m 7 — mp) 



{Tr (3Vh D V^hu) + r q ) (V ] M U V 



7/3 



4 m 2 +m 2 p 
v 2 m 2 — m 2 



7/3 



(105) 



As before, we use the third generation dominance and get the following RGEs for V a n\ 



at 



tan/? (3\V tb \ 2 h b h t + r 
4 



VM d vA V 



cd 



+ 



vM d vn v t 



at 



Id 



m t 



(VM 2 d V ] ) V cd + (VM 2 VA V td 

\ / UC v / ut 

- (3\V tb \ 2 h b h t + ; 

tanp v 



v us (y^M u v) sd ^ V uh (V^M U V 



hd 



m b 



4 r 



V us (v^M 2 u V) sd + V ub (V^M 2 V 



bd 



lQ^^-V 
at 



tan/3 (3\V tb \ 2 h b h t + r 
1 \VM 2 VA V CS + (VM 2 V^) V t 

\ / UC ' / lit 

l —h>\V tb \ 2 h b h t + ' 

n H \ 



VM d V^ V r . [VM A vA Vt 



+ 



at 



ts 



m t 



tan (3 



v ud (y^M u v) ds ^v ub (v^m u v 



bs 



m b 
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-v ud (y^M 2 u v) ds + v ub (v^m 2 u v) 



bs 



16tt 2 -K6 ^ -t&n(3(3\V tb \ 2 h b h t + r q ) 



(vM d yi) u v cb (vM d yi) u v tb 



+ 



167r 2 -V; d ~ -tanp(3\V tb \ 2 h b h 



t + r\ 



(vM d v^) c v ud (yM d v^) c V td 



+ 



m t 



(VM^) cu V ud+ (VM^) c V td 



tan/? 
4 



] - 7j (mb\ 2 h b h t + 



v cs (y ] M u v) sd V cb (V^M U V) 



m* 



+ 



bd 



m b 



v cs (y^M 2 u v) sd + v cb (v^m 2 v 



16tt 2 -V; s ^ -t&np{z\V tb \ 2 h b h 



t + r t 



(VM d vA V us (VM d V*) V t8 



+ 



m t 



(VM 2 d V^ cu V us + (VM 2 d V^) c V ts 



tan/? 
4 



— (s\V tb \ 2 h b h t 



Vcd (y ] M u v) ds v cb jy^M u v) 



+ 



6s 



m b 



-Vcd {y ] M 2 u v) ds + v c6 (v^m 2 v 



bs 



lQ7i 2 -V cb 
at 



t&n(3(3\V tb \ 2 h b h t + 



(vM d yi) cu V ub (vM d yi) ct v tb 



+ 



m. 



m t 



1 



+- 



m b tan j3 
4 



(VM 2 d V^ cu V ub+ (VM 2 d V^) c V tb 

3\Vtb\ 2 hh t + r q ) [V cd (y^M u v) db + V cs (y*M u V 



9 d 

16n 2 -V td ~ 
at 



tan /3 



K d (0M u V) d6 + \/ cs (0M M V) ; 
(3|^| 2 /i 6 ^ + r,) [(FM d l/t) ta V ud + (^M^ f ) tc V cd 
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+ 



(VMjV^) tu V ud + (VMIV^) V cd 

(3\V tb \ 2 h b h 
tanp v 



H + r t 



v ts (y^M u v) gd v tb (v^m u v 



bd 



m. 



d 
I67v 2 -V ts 
dt 



tan (3 

m t 
4 

+— 

v d 



v ts (y^M 2 u v) sd + v tb (v^m 2 v 

(3\V tb \ 2 h b h t + r q ) [{VM d V ] ) tu V us + (VM d V ] ) t V c 
(VM 2 V^) tu V us + (VM 2 d V ] ) tc V c 



tan (3 
4 r 



3\V tb \ 2 h b h t + r q 



v td (y^M u v) ds v tb (v^m u v 



bs 



m b 



lQ7r 2 -V tb 
dt 



tan p 

m t 
4 

+— 

n 



-V td (v^M 2 v) ds + V tb (V^M 2 V 
(3\V tb \ 2 h b h t + r q ) [(VM d V*) V ub + (VM d V*) V c 



(VM 2 V^ tu V ub+ (VM 2 Vl) t V cb 

} - a (mrfhbht + r q ) \V td (V^M U V)^ + V ts (v^M u V 



s6 



+- 



V td (v^M 2 u V) db + V ts (V^M 2 V 



,s6 



(106) 



We have presented the results for these RGEs even though they look quite messy because 
we believe this is the first time such an analysis has been carried out in the SUSYLR 
model, and these analytical results at the one-loop level may be useful later for future 
work in this direction. 

In order to solve these mass and mixing RGEs numerically, we need to know the 
initial values for all the 23 variables (12 masses, 9 CKM elements and 2 VEVs). We know 
the experimental values at jl = Mz for all of them except for the Dirac neutrino masses 
m Ni . We fix these values by iterations using the GUT-scale predicted values, m N .(M G ), 
which, in turn, are determined completely in terms of the other fermion masses at the 
GUT-scale in SO (10) GUT models. Here we note that adjusting the GUT-scale values 
of m^r- to fit the 5*0(10) model prediction do not change the other fermion masses at 
this scale significantly even though they are all coupled equations because of the mild 
running of the neutrino masses. Hence the mass and mixing values given in Eqs. (133]) can 
be considered as generic and independent of the specific 5*0 (10) model chosen. 

We also have the free parameters r q and r t corresponding to the couplings //* 
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and . Assuming the couplings \i a to be the same V a = 1,2,3, we have 



c=l 



1, 2, we have 



Further assuming /x* 6 = ^ V a, 6=1,2 and /i^ c = /i; V g 

r 9 = 24|/^| 2 , n = 6|^| 2 

where /i^ and can take values between and 1 (for the theory to remain perturbative). 

For the running behavior shown in Figures 2, we have chosen /i^ = 0.01 and 
Hi = 0.46 (requiring 6 — r unification) and the initial values of the Dirac neutrino masses 

m Nl (M R ) = 0.0031 GeV, m N2 (M R ) = 0.2825 GeV, m Ns = 71.86 GeV 

such that the masses evaluated at the GUT-scale, ttin^Mq), agree with those predicted 
from the specific 5*0(10) model described in Section 6. For consistency check, we note 
that the S'O(IO) model predicted eigenvalues of M D given by Eq. 

m predicted = ( _ 0028; q.2538, 77.8046) GeV, 
agree quite well with those obtained from the RGEs, 

tt4 g (M g ) = (0.0028, 0.2538, 77.8106) GeV. 
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